Half- vortex sheets and domain-wall trains of rotating two-component Bose-Einstein condensates in 
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We investigate half-vortex sheets and domain-wall trains of rotating two-component Bose-Einstein conden- 
sates in spin-dependent optical lattices. The two-component condensates undergo phase separation in the form 
of stripes arranged alternatively. The vortices of one component are aligned in lines in the low-density regions 
and filled with the other component, which results in a stable vortex configuration, straight half-vortex sheets. A 
train of novel domain walls, with spatially periodic "eyebrow-like" spin textures embedded on them, are formed 
at the interfaces of the two components. We reveal that these spatially periodic textures on the domain walls 
result from the linear gradient of the relative phase, which is induced by the alternating arrangement of the 
vortex sheets in the two components. An accurate manipulation of the textures can be realized by adjusting the 
intercomponent interaction strength, the rotating angular frequency and the period of the optical lattices. 
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I. INTRODUCTION 

In recent years, the study on topological defects has become 
a fascinating topic in Bose-Einstein condensates (BECs). In 
a single-component BEC, topological defects manifest them- 
selves as integer vortices lU-Q]. Multi-component BECs, 
which are described by a vector order parameter, allow the 
existence of more variety of exotic topological defects, such 
as fractional vortices JH-j^l, domain walls Ipl fioll and textures 
JH [TlUTill . As the simplest example of the multi-component 
condensates, two-component BECs have also attracted much 
interest to study various topological defects. 

There are various vortex configurations in rotating two- 
component BECs, such as triangular vortex lattices, rectan- 
gular vortex lattices, vortex sheets, rotating droplets and giant 
vortices lfl5U2lll . Different configurations depend on the in- 
tracomponent interaction, the intercomponent interaction and 
the particle numbers. Especially, when the intercomponent 
interaction is strong enough and the imbalance of the intra- 
component parameter is small, "serpentine" vortex sheets can 
be formed fm l20ll . However, this configuration is disordered 



and not well controlled. Although st raig ht vortex sheets were 
also discussed in the previous work ifHBoll , stable ones have 
never been obtained. This is because there exist different 
metastable vortex sheet configurations with almost the same 
energy as the straight one. 

Due to the spin degrees of freedom, the two-component 
BECs can be considered as a novel magnetic material. When 
the condensates undergo phase separation, domain walls are 
formed naturally at the interfaces of the two components. 
There have been several studies of domain walls in two- 
component BECs 122142811 . Most of them concentrate on the 
formation and dynamics of domain walls, while the internal 
structure of domain walls has not been explored so much. 
In contrast, in common magnetic materials, various internal 
structures of domain walls have been extensively investigated 
both theoretically and experimentally (see, e.g., Ref. l29ll 



and references therein). In order to study domain walls in 
two-component BECs, phase separation is required. Usually, 
phase separation is realized by the strong intercomponent re- 
pulsion. As the experimental realization of the spin-dependent 
optical lattices I3(xl3lll . the two-component condensates can 
be separated in an arbitrary form, which provides us a well- 
controlled platform to study domain walls. 

A variety of spin textures are discovered in two-component 
BECs, such as skyrmions l32tl . meron-pairs |33j] and spin-2 
textures ll34ll . The motifs of different textures always corre- 
spond to different vortex configurations. For example, an ax- 
isymmetric vortex and a nonaxisymmetric one correspond to 
a skyrmion and a meron-pair, respectively Pal , and a spin- 
2 texture corresponds to a pair of vortices with opposite sign 
that reside in different components lf34ll . This suggests that it 
is a feasible method to produce novel textures by controlling 
the arrangement of the vortices. 

In this paper, we investigate topological defects of rotat- 
ing two-component BECs in spin-dependent optical lattices. 
We find that this system supports a new stable vortex config- 
uration, straight half-vortex sheets. A train of novel domain 
walls are formed at the interfaces of the two components. 
We concentrate on their unique internal structures, and find 
that spatially periodic "eyebrow-like" spin textures are em- 
bedded on the domain walls. We reveal that these spatially 
periodic textures are directly determined by the arrangement 
of the straight half-vortex sheets. The influences of the system 
parameters on the textures are also investigated both analyti- 
cally and numerically. Our results show that the number of 
the textures on a wall is proportional to the rotating angular 
frequency and the period of the optical lattices, and with the 
increase of the intercomponent interaction strength, the tex- 
tures become thinner and thinner. This allows us to make an 
accurate manipulation of the textures. 

This paper is organized as follows. In Sec. II, we give the 
model of the rotating two-component BECs in spin-dependent 
optical lattices. In Sec. Ill, we study the stable straight half- 
vortex sheets, and discuss the discontinuity of the tangential 



component of the superfluid velocity across the sheet. In Sec. 
IV, we investigate the internal structures of the domain walls 
and reveal the formation mechanism of the spatially periodic 
"eyebrow-like" spin textures. In Sec. V, we focus on the tex- 
ture control. We conclude this paper in Sec. VI. 



II. THE MODEL 

We consider a two-level 87 Rb BEC system with 
|F = l,m/ = -l) = |1) and |F = 2,m/ = l) = |2) 
ll36ll . In the weak interaction limit, the two-component con- 
densates in a frame rotating at an angular frequency S7 around 
the z axis can be described by the coupled Gross-Pitaevskii 
(GP) equations 
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where is the macroscopic wave function of the ith com- 
ponent (i = 1,2). gij = Anh 2 aij/m represents the strength 
of interatomic interactions characterized by the intra- and in- 
tercomponent s-wave scattering lengths o,j and the mass m 
of an atom. L z = —ih(xd y — yd x ) is the z-component 
of the angular momentum operator. The external poten- 
tial consists of two parts, the harmonic trapping potential 
Vu = \m \u\ {x 2 + y 2 ) +ui 2 z 2 ~\ and the spin-dependent 
optical lattice potential Voli, where Voli = Iosin 2 (kx) and 
Vol2 = Io cos 2 (fca;). Here k is the wave vector of the laser 
light used for the optical lattice potentials and Iq is the poten- 
tial depth of the lattices. The wave functions are normalized 
as J \^i\ 2 dr = N, where N is the total number of con- 
densate atoms. 

For simplicity, we assume that the harmonic trapping fre- 
quencies satisfy uj z 3> cj±. Then, the condensates are pressed 
into a "pancake". This allows us to reduce Eq. (Q~|i to a two- 
dimensional form as [37] 
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where r\ = (h/moj z ) 1 is a reductive parameter. 
The two-dimensional wave functions are normalized as 
Si / dxdy = N. The harmonic trapping potential is 
reduced to its 2D form Vu = \miij\ (x 2 + y 2 ), where the 
tilde will be omitted in the following for simplicity. 

In order to describe the system more clearly, the intuitive 
pictures of the external potentials are presented in Figs. 1(a)- 
1(f). Experimentally, the spin-dependent optical lattice po- 
tentials Voli and Vol2 can be realized by employing two 
counter-propagating blue-detuned laser beams with the same 
frequency but perpendicular linear polarization vectors l3(ill . 
A schematic of the spin-dependent optical lattices is presented 
in Fig. 1(g). 
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FIG. 1: (Color online) (a) The harmonic trapping potential Vh is 
rotated about the z axis at a frequency fl. (b) The spin-dependent 
optical lattice potential Voli = Io sin 2 (kx), which is experienced 
by the |1) state, (c) The spin-dependent optical lattice potential 
Vol2 = Io cos 2 (kx), which is experienced by the |2) state, (d) 
The composite potential Vh + Voli- (e) The composite potential 
Vh + Vol2- (f) Cross sections of Vh (solid line), Vh + Voli (dashed 
line) and Vh + Vol2 (dotted line) along the x axis, (g) Schematic 
of the spin-dependent optical lattices. Two polarized standing wave 
laser fields a + (red) and a~ (blue) are formed by two counter- 
propagating blue-detuned laser beams with the same frequency but 
perpendicular linear polarization vectors. This gives rise to the op- 
tical lattice potentials Voli and Vol2, which are experienced by the 
|1) state (pink) and the |2) state (green), respectively. 



In our simulations, the 87 Rb atoms are assigned to the two 
states equally and the total number of them is /V = 10 5 . The 
radial and axial trapping frequencies are uj± = 2ir x 15 Hz and 
uj z = 2ir x 150 Hz, respectively. We use the scattering lengths 
lf36h : an = 100.40ao,a 2 2 = 95.00a and a i2 = a 2 \ = 
97.66eto (ao is the Bohr radius), except when we discuss the 
influence of the intercomponent interaction on the textures in 
Sec. V. The intensity of the laser light used for the optical 
lattice potentials is chosen as Io = ks x 27.5 nK {ks is the 
Boltzmann's constant) l38tl . which is powerful enough such 
that the two states are phase separated. 



HI. HALF-VORTEX SHEETS 

In this section, we present a stable vortex configuration, 
straight half-vortex sheets. By using the imaginary-time prop- 
agation method l39ll . we solve Eq. (fSJ) numerically and obtain 
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FIG. 2: (Color online) The ground state density profiles of l^il 2 , 
|i/i2| 2 and \i[)i\ 2 + \ip2 1 2 for the rotating angular frequency SI = 
0.6ix>x and the period of the optical lattice potential T = 7r£ with 
^ = (h/rrujj±) 1 ^ 2 . The locations of the vortices are marked by 
crosses (x). 



the ground state of the two-component condensates. The den- 
sity profiles with the rotating angular frequency SI = Q.6uj± 
and the period of the optical lattice potential T = 7r£ are pre- 
sented in Fig. 2, where £ = (h/muj±)~~ 1 is the spatial scale. 
All the vortices are denoted by crosses ( x ), whose positions 
are determined by the singularities of the phase. From Fig. 
2, we can see that the two-component condensates undergo 
phase separation in the form of stripes arranged alternatively. 
The vortices of one component are aligned in lines in the low- 
density regions and filled with the other component. This re- 
sults in alternatively arranged straight vortex sheets in the two 
components. Obviously, all the positions of the vortices in one 
component are vortex-free regions in the other component, so 
all the vortices are half quantized ||7H9[]. We refer to this vor- 
tex configuration as straight half- vortex sheets. Even though 
the half-vortex sheets do not influence the total density distri- 
bution of the condensates, they are crucial for the formation 
of the spatially periodic spin textures on the domain walls. 

In the absence of the spin-dependent optical lattices, the 
straight vortex sheets configuration has also been discussed in 
the phase-separated region ifTil l20ll . However, this configu- 
ration is ustable in that case, and there exist many different 
shapes of metastable vortex sheet configurations with almost 
the same energy as the straight one. This is because the energy 
of the vortex sheets is mainly determined by the intervortex 
spacing within a vortex sheet and the intersheet spacing rather 
than the shape of the vortex sheets. In contrast, when the spin- 
dependent optical lattices are present, the shape of the vortex 
sheets becomes an important factor in determining the energy. 
Any bending in the straight vortex sheets will cost much en- 
ergy, so the straight half-vortex sheets configuration obtained 
in our system is stable. In our imaginary-time propagations, 
for any trial initial configurations, the straight vortex sheets 
configuration is always uniquely obtained after sufficient con- 
vergence of the energy. 

The most well-known character of a vortex sheet is the dis- 
continuity of the tangential component of the velocity across 
the sheet. The regular arrangement of the straight vortex 
sheets allows us to observe this phenomenon clearly. In Fig. 
3(a), the tangential components (the components along the y 
axis) vi y and v^ y of the superfluid velocities vi and V2 of 
the two states are presented. We can see that both V\ y and V2 V 
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FIG. 3: (Color online) (a) The tangential components v\ v (left) and 
V2 V (right) of the superfluid velocities vi and V2 of the two states 
for SI = 0.6ojx and T = 7r£. (b) Section views of vi y (solid line) 
and viy (dashed line) along the x axis. The y component v'y=Qx of 
the rigid body rotation velocity v lb is shown by the dotted line for 
comparison. 



discontinuously jump across every sheet. In order to describe 
the tangential velocities in detail, the section views of v\ y and 
V2 y along the x axis are shown in Fig. 3(b). The y component 
a* = fix of the rigid body rotation velocity v rb = fl x r 
is presented for comparison. We find that both the tangential 
velocities vi y and V2 V have a sawtoothlike change following 
the rigid-body value v*. The value of Vi y jumps 2S16 across 
the sheet in each component and then decreases fib linearly 
in a intersheet spacing b. Here, the intersheet spacing b is 
defined as the length between two neighboring sheets in the 
same component. 

The numerical results obtained above can be understood an- 
alytically. Considering that the y component of the rigid body 
rotation velocity is v* = Vlx, we suppose that the tangen- 
tial components of the superfluid velocity on both sides of the 
sheet v~ y and vf y are independent of y. According to Onsager- 
Feynman quantization condition 



2nh 

v s • dl = N v , 



(3) 



if we choose the two sides of the sheet as the integration path, 
we can obtain that the tangential velocity jump across a sheet 
is 



Av iV = V ty - V iy = 



2ttH 1 

m d v ' 



(4) 



where d v is the intervortex spacing within a vortex sheet. This 
implies that the tangential velocity jump across the sheet is 
only determined by the intervortex spacing within the sheet. 
As the mean vortex density of each component can be esti- 
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FIG. 4: (Color online) The vectorial representations of the pseudospin S projected onto the x-y plane for (a) O = 0, T = irt; and (b) 
SI = 0.6cj±, T = 7r£. The colors ranging from blue to red describe the values of S z from —1 to 1. The locations of the vortices are marked 
by black dots (•). 



mated as 



ni = ni 



hir ' 



(5) 



we can obtain that the intervortex spacing within a vortex 
sheet is 



d ^ 
" milT 

Form Eq. (0]i and Eq. ((6]), we have 



Avs. 



2QT. 



(6) 



(7) 



As the intersheet spacing b is just equal to the period of the 
optical lattice potential T, the tangential velocity jump across 
the sheet can also be expressed as 



Av iy = 2Q.b. 



(8) 



Meanwhile, in order to follow the rigid-body value v*, Vi y 
must decrease fib in a intersheet spacing. These analytical 
results agree well with the numerical simulations above. 
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FIG. 5: (Color online) The pseudospin densities S x (left), S y (mid- 
dle), and S z (right). The upper panels show the case of f2 = and 
T — 7r£, and the lower panels show the case of Q — 0.6co± and 
T — 7r£. The locations of the vortices are marked by crosses (x). 



IV. DOMAIN- WALL TRAINS 

The spinor order parameter of the two-component BECs 
allows us to analyze this system as a pseudospin- 1/2 BEC and 
take it as a magnetic system 13511 . Introducing a normalized 
complex-valued spinor x> we represent the two-component 
wave functions as ipi = y/ pt (r)Xi (r), where pr (r) is the 
total density and the spinor satisfies |xi| 2 + |X2 1 2 = 1- In 
pseudospin representation, the pseudospin density is defined 
as S = X T crx, where a is the Pauli matrix. Then we have 

s x = 2\ X i\\x2\cos{e 1 -e 2 ), (9) 

S y = -2 |xi||xa| siii (0i-0 2 ), (10) 
s z = \xi\ 2 - \X2\ 2 , (ID 



where 0j is the phase of the wave function ipi. 

As the presence of the spin-dependent optical lattices, a 
train of domain walls are formed naturally at the interfaces 
of the two components. By the pseudospin representation, we 
investigate the response of the domain walls to rotation. In or- 
der to reveal the essential influence of the rotation on the struc- 
ture of the domain walls, the non-rotating and rotating ground 
states of the two-component condensates are calculated un- 
der the same parameters. As the rotation can create an effec- 
tive harmonic centrifugal potential with frequency £ 7 l40ll, we 
change the radial harmonic trapping frequency to \J — fi 2 
in the absence of rotation. Thus, the non-rotating ground state 
density distribution is nearly the same as the rotating ground 
state density distribution, except that no vortex is created in 
the low-density regions of each component in the non-rotating 
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FIG. 6: (Color online) (a) An amplification of Fig. 4(b) in one period. FIG. 7: (Color online) The relative phase 9i — 62 for fl — 0.6lj± 
(b) Bloch sphere of the pseudospin density vector S. Values of S z and T = 7r£, where 61 and 6 2 are the phases of the wave functions 
are represented by linear levels from blue to red (—1 to 1). ^1 and ip2, respectively. 



ground state. 

The vectorial representation of the pseudospin S for the 
non-rotating and rotating ground states are presented in Figs. 
4(a) and 4(b), respectively. Correspondingly, the pseudospin 
densities S x , S y and S z are presented in Fig. 5. From Figs. 4 
and 5, we find that a train of domain walls are formed at the 
interfaces of the spin up (S z = 1) and spin down (S z = — 1) 
regions. In the absence of rotation, the magnetic moments 
on the domain walls reverse only along the x axis [see Fig. 
4(a)], and the pseudospin density S y = [see the upper pan- 
els of Fig. 5]. Therefore, these domain walls are classical 
Neel walls. In contrast, in the presence of rotation, the mag- 
netic moments on the domain walls twist and form spatially 
periodic "eyebrow-like" spin textures [see Fig. 4(b) and the 
lower panels of Fig. 5] . This results in a train of novel domain 
walls with spatially periodic textures embedded on them. 

It should be indicated that the classical Neel wall does not 
carry topological charges. However, the twist of the mag- 
netic moments makes the novel domain wall carry topological 
charges. According to the topological charge density 

q(r) = -^-e ij S-d i Sxd j S ) (12) 

we can calculate that each period of the textures on the do- 
main wall just carries one unit topological charge. As shown 
in Figs. 4 and 5, each pair of adjacent vortices from the two 
components just induce half texture. Since a vortex in one 
component is adjacent to two vortices in the other component, 
on average, one vortex indeed induces one half unit of topo- 
logical charge. Therefore, the number of the vortices is just 
half the number of the topological charges. 

Next, we fasten our attention on the structure of the 
"eyebrow-like" spin textures. An amplification of Fig. 4(b) 
for one period is presented in Fig. 6(a). From Eq. (O and Eq. 
( [Tol l, the direction that the magnetic moments reverses along 



just depends on the relative phase, and can be represented by 
an azimuthal angle 

a = 2 (13) 

From Fig. 6(a), we can see that the azimuthal angle of the 
magnetic moments changes from — ir to tt along the domain 
wall in a period, but it is constant along the normal direction 
of the domain wall. 

It is instructive to project the pseudospin density vector 
S onto the surface of a unit Bloch sphere [see Fig. 6(b)]. 
Topologically, the topological charge counts the times that 
the Bloch sphere are covered. So one texture just covers the 
Bloch sphere once. From Fig. 6, we find that walking through 
the domain wall from one side to the other corresponds to 
strolling along a longitude line of the Bloch sphere from one 
pole to the other, and walking along the domain wall corre- 
sponds to strolling along a latitude line of the Bloch sphere. 
This is different from the case of a skymion, for which walk- 
ing along the radial direction of the skymion corresponds to 
strolling along the longitude line of the Bloch sphere, and 
walking along the azimuthal direction of the skymion cor- 
responds to strolling along the latitude line of the Bloch 
sphere. Therefore, this new texture essentially corresponds 
to a skyrmion in the polar coordinates instead of the Cartesian 
ones. 

We reveal the formation mechanism of the spatially peri- 
odic "eyebrow-like" spin textures on the domain walls. In the 
absence of rotation, there is no relative phase between the two 
components and the azimuthal angle a = 0. The magnetic 
moments on the domain walls only reverse along the x axis. 
Therefore, the domain walls are classical Neel walls. In the 
presence of rotation, straight vortex sheets are created in the 
two components and arranged alternatively on two sides of the 
domain walls. These alternatively arranged vortex sheets in- 
duce a linear gradient of the relative phase along the domain 
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FIG. 8: (Color online) The pseudospin density S y for different rotat- 
ing angular frequencies SI and periods of the optical lattice potential 
T. The upper panels show S y for T = 7r£ with (a) SI = 0Auu±, 
(b) O = 0.6ux and (c) Q, — 0.8lo±. The lower panels show S y for 
SI = 0.6cu± with (d) T = tt£, (e) T = 4/3vr£ and (f) T = 2tt£. 



FIG. 9: (Color online) (a) The pseudospin density S y for the inter- 
component scattering length 012 = 48.83an (left), 97.66an (mid- 
dle), and 976. 6ao (right) with the intracomponent scattering lengths 
an = lOCUOan and a 22 = 95.00o , SI = 0.6uj± and T = 2tt£. (b) 
Cross section views of S z along the x axis for 012 = 48.83an (solid 
line), 97.66ao (dashed line), and 976. 6an (dotted line). 



walls [see Fig. 7]. So the azimuthal angle a can be expressed 
as 



a = V (ny) , 



(14) 



where V projects the angle ny onto (— 7r, tt] and k is a constant 
coefficient, which describes the spatial change frequency of 
the azimuthal angle. The value of k will be given in the next 
section. From Eq. (fl~4l i. the azimuthal angle of the magnetic 
moments on the domain walls changes periodically along the 
y direction, and spatially periodic "eyebrow-like" spin tex- 
tures are formed. This suggests that the spatially periodic 
"eyebrow-like" spin textures on the domain walls result from 
the linear gradient of the relative phase, which is induced by 
the alternating arrangement of the straight vortex sheets in the 
two components. 



V. TEXTURE CONTROL 

In this section, we discuss how to control the "eyebrow- 
like" textures on the domain walls. Firstly, we study the dis- 
tribution of the textures analytically. The topological charge 
density q has another formulation derived from the effective 
velocity 0] 



q(r) 



2irh 



(V x v eff ) z , 



where the effective velocity is defined as 



V e ff 



(jOlVl + P2V2) 
PT 



(15) 



(16) 



with pi the density of each component. Approximately treat- 
ing the effective velocity v e ff as the classical rigid body value 
Vrb = O x r, we obtain the mean topological charge density 
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(17) 



From Eq. ( fTTI i. we can calculate that the topological charge 
on a domain wall per unit length is 



'/</ 



2nh 



(18) 



This implies that the number of the textures carried by a do- 
main wall is proportional to the rotating angular frequency 
and the period of the optical lattice potential T. From Eq. 
( fT8b . the spatial change frequency k of the azimuthal angle a 
in Eq. ( fT4l i can be calculated as 



K = 2TTl] q 



mOT 



Thus, we obtain the azimuthal angle 



a = V 



(19) 



(20) 



In order to verify the above analytical discussion, we per- 
form numerical simulations. The pseudospin density S y for 
different rotating angular frequencies SI with constant period 
of the optical lattice potential T is shown in the upper panels 
of Fig. 8, and S y for different T with constant fl is shown 
in the lower panels. Obviously, the number of the topological 
charges carried by a domain wall increases in direct propor- 
tion with the increase of SI and T. For quantitative compari- 
son, we choose O = 0.6us± and T = 7r£ as an example. From 
Eq. ( fT8l . we can calculate that the number of the textures on 
a domain wall in the region of y = [— 5£, 5£] is 3. This agrees 
well with the result of the numerical simulation in Fig. 8(b). 

In our system, as the presence of the spin-dependent opti- 
cal lattices, the two components are always phase separated 
and not subject to the immiscible condition, g\ 2 > 911(722 
l42tl . However, the intercomponent interaction g i2 plays an 
important role on the domain wall width. As the textures are 
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always concentrated on the walls, we can control the width 
of the textures by adjusting the intercomponent interaction 
gi2- The pseudospin densities S y for different intercomponent 
scattering length ai2 with constant intracomponent scattering 
lengths an and 022 are shown in Fig. 9(a). The cross section 
view of S z along the x axis is shown in Fig. 9(b). From Fig. 9, 
we can see that with increasing the strength of intercomponent 
interaction, the domain walls become narrower and narrower 
and the textures on the walls become thinner and thinner. 



VI. CONCLUSIONS 

We have investigated half-vortex sheets and domain-wall 
trains of rotating two-component BECs in spin-dependent op- 
tical lattices. A stable vortex configuration named straight 
half-vortex sheets is obtained. The discontinuity of the tan- 
gential component of the superfluid velocity across the sheet is 
discussed both numerically and analytically. We have also in- 
vestigated the response of the domain walls to rotation. In the 
absence of rotation, the domain walls are classical Neel walls 
with the magnetic moments only reversing perpendicular to 
the walls. In response to rotation, the magnetic moments on 
the domain walls twist and form spatially periodic "eyebrow- 
like" spin textures. We have revealed that these spatially pe- 
riodic textures are directly determined by the arrangement of 



the straight half-vortex sheets. The number of the textures 
carried by the domain walls can be accurately controlled by 
adjusting the rotating angular frequency and the period of the 
optical lattice, and the width of the textures can be controlled 
by adjusting the strength of intercomponent interaction. This 
allow us to make an accurate manipulation of the "eyebrow- 
like" textures. 

With the development of the magnetization-sensitive phase- 
contrast imaging technique l43ll . both the longitudinal and 
transverse magnetization of the domain and domain wall in 
BECs can be imaged non-destructively with high spatial res- 
olution fiol l44ll . We expect that the novel domain walls with 
spatially periodic "eyebrow-like" spin textures would be ob- 
served in the future experiments. 
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